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Abstract 

Theoretical and numerical studies of the Wilson flow in lattice QCD suggest that the gauge 
field obtained at flow time t > is a smooth renormalized field. The expectation values of 
local gauge-invariant expressions in this field are thus well-defined physical quantities that 
probe the theory at length scales on the order of \fi. Moreover, by transforming the QCD 
functional integral to an integral over the gauge field at a specified flow time, the emergence 
of the topological (instanton) sectors in the continuum limit becomes transparent and is 
seen to be caused by a dynamical effect that rapidly separates the sectors when the lattice 
spacing is reduced from 0.1 fm to smaller values. 



1. Introduction 

Flows in field space are an interesting tool that may allow new insights to be gained 
into the physical mechanisms described by highly non-linear quantum field theories 
such as QCD. The flow B^{t, x) of SU(3) gauge fields studied in this paper is defined 
by the equations 

Bfj, = D„G Ufl , B^\ t=Q = A^, (1.1) 

G>„ = d^B v - dvB^ + [B^ B v ], Dp = c^t + [Bp, ■ ], (1.2) 

where A^ is the fundamental gauge field in QCD (see Appendix A for unexplained 
notation; differentiation with respect to the "flow time" t is abbreviated by a dot). 
Evidently, for increasing t and as long as no singularities develop, the flow equa- 
tion (1.1) drives the gauge field along the direction of steepest descent towards the 
stationary points of the Yang-Mills action. 
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In lattice QCD, the simplest choice of the action of the gauge field U (x, fj) is the 
Wilson action [1] 



SAU) = ^]TRetr{l-t%)} 



(1.3) 



where go is the bare coupling, p runs over all oriented plaquettes on the lattice and 
U (p) denotes the product of the link variables around p. The associated flow V t (x, fi) 
of lattice gauge fields (the "Wilson flow" ) is defined by the equations 



in which d x ^ stands for the natural su(3)-valued differential operator with respect to 
the link variable Vt(x,/x) (see Appendix A). The existence, uniqueness and smooth- 
ness of the Wilson flow at all positive and negative times t is rigorously guaranteed 
on a finite lattice [2]. Moreover, from eq. (1.4) one immediately concludes that the 
action S w (V t ) is a monotonically decreasing function of t. The flow therefore tends 
to have a smoothing effect on the field and it is, in fact, generated by infinitesimal 
"stout link smearing" steps [3]. 

The Wilson flow previously appeared in ref. [4] in the context of trivializing maps 
of field space. Familiarity with this paper is not assumed, but some mathematical 
results obtained there will be used here again. An important goal in the following is 
to find out whether the expectation values of local observables constructed from the 
gauge field at positive flow time can be expected to have a well-defined continuum 
limit. Evidence for the existence of the limit is provided by performing a sample 
calculation to one- loop order of perturbation theory directly in the continuum theory, 
using dimensional regularization, and through a numerical study of the SU(3) gauge 
theory at three values of the lattice spacing. Two applications of the Wilson flow are 
then discussed, one concerning the scale-setting in lattice QCD and the other the 
question of how exactly the topological (instanton) sectors emerge when the lattice 
spacing goes to zero. 



2. Properties of the Wilson flow at small coupling 

The aim in this section partly is to show that the Wilson flow can be studied straight- 
forwardly in perturbation theory and partly to check that the expectation values of 



V t (x,fi) = -g^ {d x ^S w (V t )}V t (x, fi), 



Vt(x,fi)\ t=0 = U(x,n), 



(1.4) 
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local gauge-invariant observables calculated at positive flow time are renormalized 
quantities. 

For simplicity the perturbation expansion is discussed in the continuum theory 
using dimensional regularization. The gauge group is taken to be SU(iV) and it is 
assumed that there are Nf flavours of massless quarks. As a representative case, the 
observable 

E = \G a ^G a ^ u (2.1) 

is considered and its expectation value is worked out to next-to-leading order in the 
gauge coupling. 

2. 1 Gauge fixing 

The flow equation (1.1) is invariant under t- independent gauge transformations. E is 
therefore a gauge-invariant function of the fundamental field and its expectation 
value can consequently be calculated in any gauge. 

In perturbation theory, the gauge invariance of the flow equation leads to some 
technical complications that are better avoided by considering the modified equation 

B il = D v G vlx + \D il d v B v . (2.2) 

For any given value of the gauge parameter A, the solution of eq. (2.2) is related to 
the one at A = through 

S M = A J B M | A=0 A- 1 + Aa M A- 1 , (2.3) 

where the gauge transformation A(t,x) is determined by 

A = -Xd v B v A, A\ t=0 = 1. (2.4) 



The expectation value of E can thus be computed using the modified flow equation. 
Moreover, one is free to set A = 1, which turns out to be a particularly convenient 
choice. 

Note that the use of the modified flow equation does not interfere with the fixing 
of the gauge of the fundamental field, since E is unchanged and therefore remains a 
gauge-invariant function of the latter. 
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2.2 Solution of the modified flow equation 

In perturbation theory the gauge potential is scaled by the bare coupling, 

A» -> goAf,, (2.5) 

and the functional integral is then expanded in powers of go- The flow -B M (t, x) thus 
becomes a function of the coupling with an asymptotic expansion of the form 

oo 

£ju = ^5o-^i,fc> B tijk \ t=0 = dkiA^. (2.6) 

fc=i 

When this series is inserted in eq. (2.2), and if one sets A = 1, a tower of equations 
-B/i,fc — dvdvB^k = R^k, A; = 1,2,..., (2.7) 

is obtained, where the expressions on the right are given by 

R»,x = 0, (2.8) 
R„,2 = 2[B Vtl ,d v B li , 1 ] - [B^d^j], (2.9) 

Rn,3 = 2[B^2,^i/-B/i,i] +2[B Vi i,d v B IJli 2] 

- [B„, 2 ,d^,i] " [B^udpB^] + [B^lB^B^]}, (2.10) 
and so on. In particular, in D dimensions the leading-order equation implies 

B ll ,i(t,x) = J d D yK t (x-y)A^y), (2.11) 
f d D v ■ . 2 e" z2 / 4t 

K >^ = JwW e '" e = WF* (2 - 12) 

which shows explicitly that the flow is a smoothing operation. More precisely, the 
gauge potential is averaged over a spherical range in space whose mean-square radius 
in four dimensions is equal to V&i. 

The higher-order equations (2.7) can be solved one after another by noting that 

B fltk (t,x) = J ds J d D yK t . s (x-y)R^ k {s,y). (2.13) 



Recalling eqs. (2. 9), (2. 10), it is clear that this formula generates tree-like expressions, 
where the fundamental field is attached to the endpoints of the trees. 

2.3 Expansion of (E) 

When the series (2.6) is inserted in 

(E) = l(d„BZd„BZ-d„BZd u B«) 

+ f abc (d^B b ^B c v ) + \f abe f cde (B^BlB;Bt), (2.14) 
a sequence of terms of increasing order in g is obtained. The lowest-order term is 

So = y^B^B^ - dpB'&BU (2.15) 
and the terms at the next order are 

£i=glr bc (d,B^B b A B^), (2.16) 

£2 = g^B^B^ - 3 M B- 2 0„B- !>. (2.17) 

Each of these terms is a power series in the gauge coupling, which may be worked out 
by expressing the coefficients B^^t, x) through the fundamental field A^(x) and by 
expanding the correlation functions of the latter using the standard Feynman rules. 

In practice it is advantageous to pass to momentum space by inserting the Fourier 
representations 

B^(t,x)= [ e^e-^I», (2.18) 

J p 

B a ^ 2 (t,x) = ir bc fds f 

JO Jq,' 



e i(q+r)x e -s{q 2 +r 2 )-{t-s){q+r) 2 



X 



{<W CT - <Wa + & x (q - r)^Al{q)A- c x {r), (2.19) 



and the corresponding expressions for the higher-order fields B^^it^x). The short- 
hand 



has been introduced here in order to simplify the notation. 

Note that the flow-time integral in eq. (2.19) is similar to a Feynman parameter 
integral. In particular, if Feynman parameters are used for the diagrams contributing 
to the gluon correlation functions, one ends up with gaussian momentum integrals 
that can be easily evaluated in any dimension D. The integrals over the parameters 
then look very much the same as the ones usually encountered, except for the fact 
that the flow-time parameters are integrated up to t rather than infinity. 

2.4 Computation of £0 

In the case of the lowest-order term (2.15), the steps sketched in the previous sub- 
section lead to the formula 

So = \gl{N 2 - 1) f e- 2tp \p 2 5^-p^)D{p)^, (2.21) 

where D(p) Mil/ denotes the unrenormalized full gluon propagator. Setting D = 4 — 2e 
and choosing the Feynman gauge, the propagator assumes the form 

D (P)^ = (CP 2 V - JW^Xl - "(p))' 1 + P„P„} , (2.22) 

OO 

"(p) = X>oV) _fce "*, (2.23) 
fc=i 

from which one infers that 

«■ = ^|^> - 1) {1 + * 2( >'t^ + } • < 2 - 24 ' 

To this order, the computation is then easily completed by quoting the known result 

"1 = Tih iA7Te ~" r (& + T) - * (s - T ) + ( 2 - 25 ) 
for the gluon self-energy, 7 E = 0.577. . . being Euler's constant. 

2. 5 Computation of (E) to order 

At the next-to-leading order, the expectation value of E receives contributions from 
£ and the terms £ \ , £ 2 , . . . up to order g$ generated by the expansion of the flow 
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equation (cf. subsect. 2.3). Some of the latter involve the 3-point gluon vertex, but 
in all cases only the leading-order gluon correlation functions are required. 

In the case of the term £2, for example, the purely algebraic part of the calculation 
leads to the integral 

ft n (2t-s)p 2 -s(q 2 +r 2 ) 

£2 = N(N 2 - l)gt / ds / 

JO Jq.r p q r 

x{(D- l)p 2 (p 2 + q 2 + r 2 ) + 2{D - 2)(q 2 r 2 - (qr) 2 )} + O(^) (2.26) 

where p = q+r. This integral appears to be quite complicated, but the polynomial in 
the numerator of the integrand allows the expression to be simplified and eventually 
to be evaluated analytically (see Appendix B for further details). The result of the 
computation, 

£2 = N(N 2 - l) ^jL—- {*. - § +451n2 - f ln3 + 0(6)} + 0( 9 °), 

(2.27) 

shows that these contributions are not free of ultra-violet singularities. 

The computation of the other terms follows the same pattern and does not present 
any additional difficulties. Collecting all contributions, the result 

N 2 - 1 

<£> = boj^joj2-(D -!){! + c l9 2 + O( 5o 4 )}, (2.28) 
Cl = JL( 4 ,r)W{tf (M + f - 31n3 ) " N * (& + I ~ I ln2 ) + 

(2.29) 

is then obtained. 
2.6 Renormalization 

The bare coupling g$ is related to the renormalized coupling g in the MS scheme 
and the associated normalization mass fi according to [5] 

g 2 = <?V e (47re-^r {l - \b G g 2 + 0(/)}, (2.30) 
= JL{£tf-§JV f }. (2.31) 



Now if eq. (2.28) is written as an expansion in the renormalized coupling, the terms 
proportional to 1/e cancel and one obtains 

^ ) = 3 m^^ 1 + c " lg2 + Q(g4) ^ (2 ' 32) 

~ C1 = JL { AT (£L + f - 3 In 3) - N t (§L + § - | In 2) } , (2.33) 

at e = 0, where L = ln(8/z i) + 7e- To this order in the gauge coupling, (E) thus 
turns out to be a quantity that does not need to be renormalized and which therefore 
encodes some physical property of the theory. 

In terms of the running coupling a(q) at scale q = (8t) -1 / 2 , the expansion (2.32) 
assumes the form 

{E) = ^l^ w + ha{q) + 0(a2)} ' (2 - 34) 

fc i = J^{ N (t7e + f - 3 In 3) - N t (| 7 e + | - | In 2)}. (2.35) 
In particular, for N = 3 one obtains 

(E) = —,a(q){l + k 1 a(q) + 0(q 2 )}, fci = 1.0978 + 0.0075 x N t . (2.36) 

The next-to-leading order correction is reasonably small in this case and corresponds 
to a change in the momentum q by a factor 2 or so if Nf < 3. 



3. Lattice studies of the Wilson flow 



In QCD the perturbation expansion (2.36) is expected to be applicable at small flow 
times only, where the smoothing range y/St is at most 0.3 fm or so. The properties 
of the Wilson flow at larger values of t can however be studied straightforwardly 
using the lattice formulation of the theory and numerical simulations. 

The simulations of the SU(3) gauge theory reported in this section mainly serve 
to clarify whether (E) scales to the continuum limit as suggested by perturbation 
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Table 1. Lattice parameters, statistics and reference flow time 



Lattice 


P 


a [fm] 


Acnfg 


t /a 2 


48 x 24 3 


5.96 


0.0999(4) 


100 


2.8449(64) 


64 x 32 3 


6.17 


0.0710(3) 


100 


5.579(14) 


96 x 48 3 


6.42 


0.0498(3) 


100 


11.364(24) 



theory. Along the way, a new scale-setting method will be proposed based on the 
observed properties of (E). 

3.1 Simulation parameters 

Ensembles of representative gauge fields were generated on three lattices, using the 
Wilson gauge action and a combination of the well-known link-update algorithms. 
The values of the lattice spacing quoted in table 1 derive from the results in lattice 
units for the Sommer reference scale r = 0.5 fm [6] published by Guagnelli et al. [7]. 
At the chosen couplings j3 = 6/<?q, the spacings of the three lattices thus decrease 
from roughly 0.1 to 0.05 fm by factors of l/\/2. Moreover, the lattice sizes in physical 
units are approximately constant. 

In order to safely suppress any residual statistical correlations of the iV cn f g gener- 
ated fields, the separation in simulation time of the fields was taken to be at least 10 
times the integrated autocorrelation time of the topological charge. The definition 
of the latter on the lattice is ambiguous to some extent, but its autocorrelation time 
is largely independent of the choices one makes and is known to increase very rapidly 
when the lattice spacing is reduced [8,9]. In particular, already at a = 0.07 fm all 
other usual quantities of interest tend to be far less correlated in simulation time. 

3.2 Observables 

For any given gauge field configuration U (x, [/,), the flow equation (1.4) can be inte- 
grated numerically up to the desired flow time and one may then construct gauge- 
invariant local observables from the gauge field at this time. In particular, 

E = 2 ^ Reti{l - V t (p)} (3.1) 

p€P x 

is a possible definition of the density E on the lattice, P x being the set of unoriented 
plaquettes with lower- left corner x. 
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Fig. 1. The anti-hermitian traceless part of the average of the four plaquette Wilson 
loops shown in this figure can be taken as the definition of the field tensor a 2 G, lv {x) 
on the lattice. All loops are in the (/x, i/)-plane, have the same orientation and start 
and end at the point x. 



Another more symmetric definition of E is obtained by introducing a lattice ver- 
sion of the field tensor G^ v {x) (see fig. 1). E is then simply given by the continuum 
formula (2.1). Both definitions are equally acceptable at this point, since they re- 
spect all formal requirements (locality and gauge invariance in particular) and since 
they converge to the correct expression in the classical continuum limit. 

For the numerical integration of the flow equation (1.4), any of the widely known 
integration schemes can in principle be used (see ref. [10], for example). The Euler 
integrators previously discussed in ref. [4] are particularly simple but also the least 
efficient ones. Evidently, the integration errors should be much smaller than the 
statistical errors of the calculated quantities. A fairly simple and numerically stable 
integrator that allows this condition to be easily met is described in Appendix C. 

3.3 Time dependence of {E) 

To leading-order perturbation theory, the dimensionless combination t 2 (E) is a con- 
stant proportional to the gauge coupling. At the next order, the scale invariance of 
the theory is broken and t 2 (E) develops a non-trivial dependence on the flow time 
t. Asymptotic freedom actually implies that the combination slowly goes to zero in 
the limit t — > 0. 

Recalling the result 



Al^,, = 0.602(48)/r (3.2) 



for the A-parameter in the MS scheme obtained by the ALPHA collaboration [11], 
and using the four-loop evolution equation for the running coupling a{q) [12], the 
perturbation series (2.36) can be evaluated at any value of the flow time given in 
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Fig. 2. Simulation data obtained at a = 0.05 fm for t 2 (E) as a function of the flow 
time t (black line). Statistical errors are smaller than 0.3% and therefore invisible on 
the scale of the figure. The curve predicted by the perturbation series (2.36) and the 
known value (3.2) of the A-parameter is also shown (grey band). 

units of tq. The curve obtained in this way is shown in fig. 2 together with the error 
band that derives from the error of A quoted in eq. (3.2). 

Perhaps somewhat fortuitously, the simulation results obtained at a = 0.05 fm 
accurately match the perturbative curve over a significant range of t. The symmetric 
definition of E has here been used and for clarity the data from only one lattice are 
shown (as discussed below, the lattice-spacing effects are small and the data from 
the other lattices would therefore lie nearly on top of the line plotted in fig. 2). 

Beyond the perturbative regime, t 2 {E) grows roughly linearly with t, at least so 
within the range covered by the simulation data. The slowdown of the density (E) 
from the perturbative 1/t 2 to a smoother \/t behaviour may perhaps be explained 
by noting that the Wilson flow tends to drive the gauge field towards the stationary 
points of the gauge action. In the vicinity of these points of field space, the right- 
hand side of the flow equation (1.4) is small and E consequently changes only little 
with time. 

3.4 Lattice-spacing effects 

Perturbation theory suggests that the density (E) scales to the continuum limit like 
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a = 0.05 fm 0.07 fm 



0.1 fm 



i I 



0.96 - 




0.90 - 







0.01 



0.02 



0.03 



0.04 



(air,) 1 



Fig. 3. Extrapolation of the dimensionless ratio \/8to/ro of reference scales to the 
continuum limit (open data points). The black data points were obtained using the 
symmetric definition of E and the grey ones using the expression (3.1). 

a physical quantity of dimension 4. The scaling behaviour of (E) can be checked by 
introducing a reference scale to through the implicit equation 



(see fig. 2). If (E) is physical, the dimensionless ratio to/rQ must be independent of 
the lattice spacing, up to corrections vanishing proportionally to a power of a. 

The data plotted in fig. 3 clearly show that the ratio of reference scales smoothly 
extrapolates to the continuum limit. As may have been suspected, the lattice effects 
appear to be of order a 2 . They are at most a few percent on the lattices considered 
and particularly small if the symmetric definition of E is employed. At other points 
in time, the scaling violations behave similarly, but tend to increase towards small t, 
where the smoothing range y&t is only 2 or 3 times larger than the lattice spacing. 

3.5 Synthesis 

(a) Continuum limit. The numerical studies of the Wilson flow strongly support the 
conjecture that the gauge fields generated at positive flow time are smooth renor- 
malized fields (except for their gauge degrees of freedom). In the case of QCD with 
a non-zero number of sea quarks, similar studies however still need to be performed. 



(3.3) 



12 



Additional confirmation in perturbation theory and through the consideration of 
further observables is evidently desirable. 

(b) Scale setting. The time to defined through eq. (3.3) may serve as a reference scale 
similar to the Sommer radius rg. With respect to the latter, to has the advantage that 
its computation does not require any fits or extrapolations. Moreover, an ensemble 
of only 100 independent representative field configurations allows t to be obtained 
with a statistical precision of a small fraction of a percent (for illustration, the values 
of to /a 2 computed using the symmetric definition of E are listed in table 1). 

(c) Universality. On the lattice the definition of local gauge-invariant quantities like 
E is not unique, but the results obtained in this section indicate that the differences 
become irrelevant in the continuum limit (see fig. 3). This kind of universality is a 
consequence of the fact, further elucidated in sect. 4, that the gauge fields generated 
by the Wilson flow are smooth on the scale of the lattice spacing. The universality 
classes of the local fields composed from the gauge field at fixed t/to are therefore 
determined by the asymptotic behaviour of the fields in the classical continuum 
limit. 



4. Functional integral and topological sectors 

In lattice gauge theory, the space of gauge fields is connected and the concept of a 
topological sector has therefore no a priori well-defined meaning. However, one also 
knows since a long time that any classical continuous gauge field (including multi- 
instanton configurations) can be approximated arbitrarily well by lattice fields. The 
sectors are hence included in the field space but are not separated from one another. 

An understanding of how exactly the sectors get divided when the lattice spacing 
is taken to zero can now be achieved using the Wilson flow. The existence of the 
topological sectors thus turns out to be a dynamical property of the theory rather 
than being a consequence of an assumed or imposed continuity of the fields. 

4-1 Field transformation 

On a finite lattice, however large, the transformation 

U -> V = V to (4.1) 

is invertible and actually a diffeomorphism of the field space [4] . One can therefore 
perform a change of integration variables in the QCD functional integral from the 
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fundamental field to the field at flow time t . As already noted in ref. [4], the asso- 
ciated Jacobian can be worked out analytically and be expressed through the Wilson 
action along the flow. The expectation value of any observable 0(V) is then given 



where S(U) denotes the total action (including the quark determinants) of the theory 
before the transformation. 

Evidently, the functional integral may be rewritten in this way as an integral over 
the gauge field at any flow time t. Setting t to the reference time to is however an in- 
teresting and natural choice. In particular, the fields that dominate the transformed 
integral then have a characteristic wavelength on the order of the fundamental low- 
energy scales of the theory. 

4-2 Smoothness of the dominant fields 

A quantitative measure for the smoothness of a given lattice gauge field V is 

/i = maxs„, s p = Retr{l — V(p)}, (4.4) 
v 

the maximum being taken over all plaquettes p (as before, V(p) denotes the product 
of the link variables around p). In the functional integral (4.2), all possible values 
of h occur, but large values are expected to be unlikely in view of the fact that the 
expectation value (s p ) scales proportionally to a 4 . 

Since all terms in the action (4.3) have the same sign and favour smooth configura- 
tions, it is certainly plausible that large plaquette values s p are strongly suppressed. 
The probability for s p on a given plaquette p to be larger than some specified value 
s in fact decreases roughly like a 10 when the lattice spacing is reduced (see fig. 4). In 
a finite volume of fixed physical size, the statistical weight of the field configurations 
with values of h beyond a given threshold is therefore rapidly going to zero in the 
continuum limit, while the fields that dominate the transformed functional integral 
become uniformly smooth in space. 



by 




(4.2) 




(4.3) 
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s 

Fig. 4. Plot of the probability for s p [eq. (4.4)] on a given plaquette p to be above 
some specified value s. From top to bottom, the three curves shown correspond to the 
values a — 0.1, 0.07 and 0.05 fm of the lattice spacing. The calculation was performed 
in the SU(3) gauge theory using the ensembles of fields described in subsect. 3.1. 

4-3 Dynamical separation of the topological sectors 

Many years ago, the space of lattice gauge fields was shown to divide into discon- 
nected sectors if only fields satisfying a certain smoothness condition are admitted 
[13,14]. The proof is based on a geometrical construction of a local lattice expression 
for the topological charge which assumes integer values and which has the correct 
classical continuum limit. 

In the case of QCD, the fields satisfying f 

h < 0.067 (4.5) 

are included in the subspace covered by the geometrical construction and thus fall 
into topological sectors very much like the continuous gauge fields in the continuum 
theory. The space of fields "between the sectors" may accordingly be characterized 
by the inequality h > 0.067. 

f The theorems proved by Phillips and Stone [14] assume a simplicial lattice. To be able to apply 
them in the present context, the hypercubic cells of the lattice must be divided into simplices. An 
interpolation of the gauge field to the added links is then required and can easily be accomplished 
using the same interpolation methods as the ones applied elsewhere in refs. [13,14]. 
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Recalling the discussion in subsect. 4.2 of the smoothness properties of the fields 
that dominate the functional integral (4.2), the topological sectors are now seen 
to emerge as a consequence of the fact that the fields with large values of h have a 
rapidly decreasing weight in the integral when the lattice spacing is taken zero. In the 
case of the simulations reported in sect. 3, for example, the fraction of representative 
fields satisfying the bound (4.5) increases from 0% at a = 0.1 fm to 8% and 70% 
at a = 0.07 and 0.05 fm, respectively. The asymptotic behaviour of the weight of 
the fields between the sectors cannot be safely determined from these data, but the 
curves shown in fig. 4 suggest that the probability of finding such a configuration on 
lattices of a given physical size goes to zero proportionally to a 6 . 

4-4 Topological susceptibility 

Once the functional integral (4.2) splits into a sum of integrals over effectively discon- 
nected sectors, the assignment of the topological charge Q to the field configurations 
in a representative ensemble of fields becomes unambiguous. The geometrical con- 
struction mentioned before could be used to compute the charge, but in view of the 
universality property discussed at the end of sect. 3, a straightforward discretization 
of the topological density, using the symmetric lattice expression for the field tensor 
G M „, is expected to give the same results for the moments (Q n ) in the continuum 
limit. 

On the lattices listed in table 1, the second moment (Q 2 ) computed along these 
lines is about 50 and the topological susceptibility Xt turns out to be independent of 
the lattice spacing within statistical errors. A fit by a constant then gives the result 

x] /4 = 187.4(3.9) MeV, (4.6) 

which differs by less than two standard deviations from the value 194.5(2.4) MeV 
[16] obtained at flow time t = using a chiral lattice Dirac operator and the index 
theorem [15]. Moreover, as shown in a forthcoming publication, the ratios of the 
fourth root of the second moments computed here and through a variant [18] of the 
universal formula proposed in [17] extrapolate to 1 in the continuum limit to within 
a statistical uncertainty of 3%. 

All these empirical results support the conjecture that the moments of the charge 
distribution in the functional integral (4.2) coincide with those obtained from the 
index theorem [15-18] and thus the ones appearing in the chiral Ward identities 
[19,20]. However, while highly plausible, the equality remains to be theoretically 
established. 
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5. Concluding remarks 

The results reported in this paper shed some new light on the nature of non-abelian 
gauge theories and the continuum limit in lattice QCD. They are obviously incom- 
plete in several respects and give rise to some interesting questions. In particular, an 
all-order analysis of the Wilson flow in perturbation theory will probably be required 
in order to achieve a structural understanding of its renormalization properties. 

The one-loop calculation in sect. 2 suggests that the fields obtained at positive flow 
time are renormalized fields for any number of quark flavours. So far the question 
was studied numerically only in the pure gauge theory, but it may be encouraging 
to note that the conjecture is true in QED, for any charged matter multiplet, as a 
consequence of the gauge Ward identity. 

An intriguing aspect of the transformed functional integral (4.2) is the fact that it 
is dominated by smooth fields. Since the field transformation (4.1) is invertible, the 
integral nevertheless encodes all the physics described by the theory. Some properties 
(the division into topological sectors, for example) however become particularly 
transparent in this formulation of the theory, while its behaviour at high energies is 
better discussed in terms of the fundamental gauge field. 

I wish to thank Peter Weisz for a critical reading of the paper and for checking 
some of the formulae in sect. 2. All numerical simulations were performed on a ded- 
icated PC cluster at CERN. I am grateful to the CERN management for providing 
the required funds and to the CERN IT Department for technical support. 



Appendix A. Notational conventions 

The Lie algebra su(N) of SXJ(N) may be identified with the linear space of all anti- 
hermitian traceless N x N matrices. With respect to a basis T a , a = 1, . . . , A 2 — 1, 
of such matrices, the elements X 6 su(A) are given by X = X a T a with real compo- 
nents X a (repeated group indices are automatically summed over). The structure 
constants f abc in the commutator relation 

[ T a )T b ] = f abc T C (A.l) 

are real and totally anti-symmetric in the indices if the normalization condition 
tr{T a T b } = -\5 ab (A.2) 
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is imposed. Moreover, f acd f bcd = N5 ab . 

Gauge fields in the continuum theory take values in the Lie algebra of the gauge 
group. Their normalization is usually such that gauge transformations and covariant 
derivatives do not involve the gauge coupling (see, however, subsect. 2.2). Lorentz 
indices ft, u, ... in D dimensions range from to D — 1 and are automatically summed 
over when they occur in matching pairs. The space-time metric is assumed to be 
Euclidian. In particular, p 2 = p^p^ for any momentum p. 

The lattice theories considered in this paper are set up as usual on a four-dimen- 
sional hyper-cubic lattice with spacing a. In particular, lattice gauge fields U(x, fx) 
reside on the links (x, fi) of the lattice and take values in the gauge group. The link 
differential operators acting on functions f(U) of the gauge field are 

d (T a if (y,u) = (x,u), 

d a x J{U) = ^-f{e sX U) , X(y,u) = \ (A.3) 

as 8=0 1 otherwise. 

While these depend on the choice of the generators T a , the combination 

d x ,,f(U) = T a d a x J(U) (A.4) 

can be shown to be basis-independent. 



Appendix B. One-loop integrals 

The Feynman integrals other than £ which contribute to (E) at order g$ involve an 
integration over two momenta, q and r, and over at most two flow-time parameters. 
In all cases the integrand is of the form 



q 2 r 2 (q + r) 2 



(B.l) 



where P(q, r) is a Lorentz-invariant polynomial and s, u, v are linear combinations 
of the flow-time parameters. 

As already mentioned, these integrals can in principle be evaluated by substituting 
the Feynman parameter representation for the propagators 1/q 2 , 1/r 2 and l/(q+r) 2 . 
A more economic computation is, however, always possible using the basic integrals 

[- sq > ( o,- a _ s« T(D/2-a) 

J q iq) ~WW 2 T{D/2) ' (R2) 
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/ 5 = ,a ^nm 5w , \ ( su + uv + «*) 7 ■ (B.3) 

V Q 2 (Ak) d {D -2){u + vY 

The second formula, for example, allows the integral 



/ 

Jq,r 



t {q 2 +r 2 + (q+r) 2 ) i 

— = (S)W{ 8ta2 - 41113 + < t »> < B - 4 > 



to be quickly evaluated once the propagator 1/r 2 is replaced by its Feynman param- 
eter representation. 

A special case are integrals like 

* r e -{t+s)(q 2 +r 2 )-(t-s)(q+r) 2 

ds / — qr 

o Jg,r Q 2 r 2 

= (4^)W { "" 41n2 + 21n3 + 0(e)} ' (R5) 
whose integrand is proportional to qr. Noting 

qr=±{(q + r) 2 -q 2 -r 2 }, (B.6) 

the factor can be traded for a differentiation with respect to the flow-time parame- 
ters. Most of the time, this allows one integral over these parameters to be performed 
right away and thus leads to simpler integrals without factors of qr. 



Appendix C. Numerical integration of the Wilson flow 

On a finite lattice, the space Q of all gauge fields is a finite power of the gauge group 
and thus itself a Lie group. The associated Lie algebra g coincides with the linear 
space of all link fields with values in the Lie algebra of the gauge group. From this 
abstract point of view, the flow equation (1.4) is an ordinary first-order differential 
equation of the form 

V t = Z(V t )V u (C.l) 

where V t G Q and Z(V t ) G g. 
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The Runge-Kutta scheme described in this appendix obtains the solution of the 
flow equation at times t = ne, n = 1, 2, 3, . . ., recursively, starting from the initial 
configuration at t = 0. The rule for the integration from time t to t + e is 

Wo = V t , 

Wi = exp{\Z }W , 

W 2 = exp{|Z 1 -^Z }Wi ) 

y t+e = exp{|Z 2 -§Z 1 + i|Zo}^ 2 , (C.2) 
where 

Z l = eZ(^), i = 0,l,2. (C.3) 

Note that this rule is fully explicit. Moreover, since the gauge field can be overwritten 
from one equation to the next, and since Z can be overwritten by |Zi — t^Zq, 
intermediate storage space for only one of these latter fields is required. 

A straightforward calculation shows that the integration scheme (C.2) is accurate 
up to errors of order e per step. The total error of the integration up to a specified 
flow time thus scales like e 3 . Empirically one finds that the integration is numerically 
stable in the direction of positive flow time, the integration errors in the link variables 
being on the order of 10~ 6 if e = 0.01. 
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